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ABSTRACT 

Aims. To constrain the equation of state of super-nuclear density matter and probe the interior composition of the X-ray pulsar in SAX 
J1808.4-3658. In our estimation, we consider both its persistent 2.49 ms X-ray pulsations discovered by Wijnands and van der Klis 
from using the Rossi X-ray Timing Explorer, which is interpreted to come from an accreting-powered millisecond X-ray pulsar in the 
low mass X-ray binaries, and the corresponding mass-radius data analyzed of the light curves of SAX J1808.4-3685 during its 1998 
and 2005 outbursts by Leahy et al. from assuming a hot spot model where the X-rays are originated from the surface of the neutron 
star. 

Methods. The interior composition of the X-ray pulsar in SAX J1808.4-3658 is obtained from comparing the mass-radius relations, 
Keplerian motions and gravitational wave radiation instabilities of the predictions under dilferent theoretical models with the analyti- 
cally mass-radius district and rotation frequency of observational data. 

Results. We show that SAX J1808.4-3658 may include the exotic matter such as hyperon or quark matter, but we can't distinguish 
the star with hyperon matter from quark matter. 
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1. Introduction 

The properties of nuclear matter under extreme conditions 
are correlative to the in-depth experiment and theoretical re- 
search of the heavy nuclear physics and nuclear astrophysics 
intimately. Neutron stars provide us a unique spherical envi- 
ronment to study the properties of matter with low temper- 
ature and high density. Since the discovery of the first pul- 
sar by Hewish et al. (1968) and then the confirmation of 
it to be a fast rotational neutron star by Gold and Pacini 
(1968), the compositions and properties of the interiors of neu- 



tron stars have attracted much attention (Pandharipande 1971 



Glendenning 1985| 


Glendenning et al. 1992, Sahu et al. 1993 


Thorsson et al. 1994 


Alford & Reddy 2003)1. If we have made 



this significant problem clear, it would modify our understand- 
ing on the structure of the matter, strengthen the cognition about 
the formation of the universe and in all contribute most to the 
research of astrophysical physics, nuclear physics and particle 
physics. 

However, in fact due to the uncertain nuclear physics, 
the equation of state for neutron star at super-nuclear 
density is still an indeterminacy, on which the maxi- 
mum mass, radius and rotation frequency depend strongly. 
Many investigators spontaneously expect to probe the 
matter under these conditions and constrain or rule out 
some equations of state through astrophysical observations 
([Glendenning & Moszkowski 1991] ILattimer & Pra kash 2004, 
Lackey et al. 2006 1 |Lavagetto et al. 2006] IKlahn et al. 20061 



Lattimer & Praka sh 20071 [Pan & Zheng 2b07| l. Although re- 
searchers universally impose the inferred masses and radii or 
rotation frequencies for pulsars on the equations of state, their 



common treatments cannot uniquely examine and distinguish 
all the classes of compositions inside neutron stars, i.e. it could 
rule out neither traditional neutron stars nor hyperon stars or 
quark stars determinately. This problem should be carefully 
considered in other effective ways for the further estimation. 



2. The X-ray pulsar in SAX J1 808.4-3658 

In this letter, we mainly estimate the X-ray pulsar in the transient 
X-ray burster SAX J1808.4-3658. Wijnands and van der Klis 
discovered its persistent 2.49 ms X-ray pulsations from using 
the Rossi X-ray Timing Explorer, and interpreted them to come 
from an accreting-powered millisecond X-ray pulsar in the low 
mass X-ray binaries (Wijnands & van der Klis 1998 ). Recently, 
Leahy et al. have analyzed its light curves during 1998 and 2005 
outbursts to get the corresponding mass-radius data from assum- 
ing a hot spot model where the X-rays are originated from the 
surface of the neutron star ( |Leahy etsl. 2008[ ). 



3. Probing stellar matter 

In what follows, we present our approach to probe the matter 
composition of the X-ray pulsar in SAX J1808.4-3658 using the 
two observational characteristics mentioned above. 



3.1. Mass-Radius estimates 

The structure of a neutron star is determined by the local bal- 
ance between the attractive gravitational force and the pres- 
sure force of the neutron star matter. The X-ray pulsar in 
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Table 1. Summary of equations of state Al-Dl introduced in 
the text, here compositions refers to the interacting components, 
n - neutrons, p - protons, e - electrons, H - hyperons and Q - 
quarks. GC and MC represent the deconfinement phase transi- 
tion under Gibbs and Maxwell constructions respectively. RMF 
is the abbreviation for relativistic field theoretical approach in 
the mean field approximation, MIT the simple MIT bag model, 
eMIT the effective mass MIT bag model, APR the equation of 
state of Akmal, Pandharipande and Ravenhall using the varia- 
tional chain summation, pQCD perturbative QCD corrections, 
and BBG the Brueckner-Bethe-Goldstone many-body approach. 



Symbols 


Approaches 


References 


Compositions 


Al 


RMF 


Glendenning (1997) 


npe 


A2 


RMF 


Lackey et al. (2006) 


npeH 


A3 


RMF+MIT 


Pan & Zheng (2007) 


npeQ (GC) 


A4 


RMF+eMIT 


Zheng et al. (2007) 


npeQ (MC) 


Bl 


APR 


Akmal etal. (1998) 


npe 


B(2-3) 


APR+pQCD 


Alford et al. (2005) 


npeQ (GC/MC) 


CI 


BBG 


Nicotra et al. (2006) 


npe 


C(2-3) 


BBG-hMIT 


Nicotra et al. (2006) 


npeQ (GC/MC) 






Fahri& Jaffe(1984) 




C4 


BBG 


Baldo et al. (2000) 


npeH 


Dl 


eMIT 


Schertler et al. (1997) 


Q(u,d,s) 



SAX J1808.4-3658 rotates at 2.49 ms, which in fact should be 
solved through the perturbative approach of the non-rotating 
equilibrium configuration in the framework of general relativ- 
ity (Hartle 1967 Hartle & Thor ne 196"8] l. But the rotation cor- 
rections of the radii and masses are about 10% and 2-3%, which 
couldn't affect our results significantly here, the non-rotating ap- 
proximation suffices. 

According to the Tolman-Oppenheimer-Volkoff theory 
(TOV) under the consideration of the effect of general rel- 
ativity, the structures of spherically symmetric static neutron 
stars could be determined by a set of equations (ITolman 19391 
Oppenheimer & Volkofi^ 19"39ll: 



dp{r) [e{r) + p(r)] [m{r) + AuPpir)] 



dr 
dr 



r\r — 2m(r)] 



(1) 



(2) 



(G = c = 1) where p{r) and e{r) are the pressure and energy 
density of the matter at the radius r, and m{r) is the total mass 
inside the star within a sphere of given radius r: 



Jo 



m(r) — A-n I e(r')r'^dr' . 



(3) 



After the equation of state of the star e(r) - e(p(r)) is given, the 
TOV equations could be solved as an initial value problem. For 
a given equation of state, there exists a unique relationship be- 
tween the stellar mass and radius. Thus different mass-radius re- 
lations obtained by considering various equations of state could 
be used to distinguish one from another. 

As we have mentioned before, the neutron stars actually 
could have many classes. Besides traditional neutron star (Al, 



Bl, CI) that mainly consists of neutrons (Glendenning 1997 
lAkmal et al. 19981 INicotra et al. 2006] l. when hyperons become 
populated in addition to the nucleons, it corresponds to the hy- 
peron star (A2, C4) (L ackey etal. 2006[ [Baldo et al. 20001 ). In 
accordance with the predictions of quantum chromodynamics, 
deconfined quark matter may also be produced there. When 



the highly compressed matter could undergo the deconfinement 
phase transition into u, d, s quark matter, or according to ab- 
solute stable strange quark matter hypothesis, the quark-hybrid 
star (A3-4, B2-3, C2-3) (|Pan & Z heng 2007, Zheng et al. 2007| 
I Alford et al. 20051 INicotra et al. 2006 Fahri & Jaffe 19841 i and 
strange star (Dl) dSchertler et al. 1997l l may exist. In Fig. 1, we 
compare various radius-mass relations with the observation data 
of SAX J1808. 4-3658, of which the equations of state have been 
listed detailedly in Table 1 . Obviously, only these softer equa- 
tions of state could satisfy the inferred mass-radius region, and 
some of the neutron stars, hyperon stars, quark stars and strange 
stars match the observation well. But as predicted above, we 
can't tell the exact components in the core of this neutron star. 



3.2. Rotation estimates 

For a uniform rigid compact star with mass M and radius R, there 
exist some limits on the attainable rotation frequency. The most 
obvious one is the Keplerian limit dLattimer & Prakash 20041) 

(M/Mq)1/2 



vk = L042X 10" X 



(4) 



(R/lQkm)^^ 

which is also called the mass-shedding limit. When its spin fre- 
quency exceeds this limit, the matter near the equator would run 
away, and the star may be no longer stable. It is foreign to the 
material qualities of the matter inside. So if we get the frequency 
of observed pulsar and assume it to be the Keplerian limit, we 
can obtain a critical mass-radius relation to constrain the equa- 
tions of state, which has been given as the cyan district in Fig. 1, 
but it also seems to have no help in distinguishing the composi- 
tion of the super-nuclear matter even if the mass-radius method 
is considered together 

Actually, the emission of gravitational radiation following 
the excitation of non-radial oscillation modes may lead to the in- 
stability of rotating stars dChandrasekhar 19'70llAndersson 19981 
lAndersson & Kokkotas 20(JTl l, and the corresponding limiting 
rotation can be obtained via 

— + -^0. (5) 

Here tc < is the characteristic time scale for energy loss due to 
gravitational waves emission, and Ty the damping timescales due 
to shear, bulk viscosities and other rubbings, which relates to the 
viscosities of the matter inside neutron stars and their differences 



could result in diverse behaviors (Pan & Kang 2007). 

In Fig. 2, we show the limit rotation frequencies for different 
classes of neutron stars under the relativistic mean field theory 
for demonstration, which is constructed by means of the method 
brought forward by Zheng et al.(2007). Clearly, these thick solid 
lines should be the genuine upper frequencies for each class. 
After matching our theoretical predictions with the observation 
data, we find that the traditional neutron star should be excluded, 
although its Keplerian limit could reach above the 401 Hz of 
SAX J1808.4-3658. While hyperon star, quark star and strange 
star are supposed to be the best candidates. Therefore, the star 
could contain either hyperons or quarks at super-nuclear region, 
but we can't tell them from each other We should emphasize 
that, the conclusion is still true if we apply it to various classes 
of neutron stars under other approaches. 

4. Conclusions 

The composition of matter in the core of neutron stars has at- 
tracted much attention owing to its important significance. Up 
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Fig. 1. Radius-mass relations for various equations of state of 
compact stars listed in Table 1. The red and cyan districts re- 
fer to the causality limit and simply traditional Keplerian rota- 
tion limit at 2.49 ms discovered by Wijnands and van der Klis 
(1998) respectively. These green contours correspond to 80%, 
90%, 95.4% and 99.7% confidence levels of the X-ray pulsar in 
SAX J1808.4-3658 analyzed by Leahy et al. (2007). Note that 
we here use the non-rotating approximation in the calculation, 
as the radii and masses of these stars are about 10% and 2-3% 
larger due to the rotation correction at 2.49 ms, which couldn't 
affect our results significantly. 



to now, SAX J 1808.4-3685 is one of seven known accreting ms 
X-ray pulsar Its persistent 2.49 ms X-ray pulsations was discov- 
ered by Wijnands and van der Klis and interpreted as due to the 
ms rotation of the center neutron star And Leahy et al. have an- 
alyzed its corresponding mass-radius relation. In our estimation, 
we try to probe the inner components of it in our own way by 
comparing the mass-radius relations and genuine rotation fre- 
quencies under different theoretical models with the observa- 
tional data. We finally come to a conclusion that the pulsar in 
SAX J1808. 4-3685 is a star containing exotic matter, which just 
bases on these two observational properties, but we can't distin- 
guish it with hyperon matter from quark matter at the present 
time, which must depend on more observational information 
about it, such as the thermal emission data. 
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Fig. 2. Limit rotation frequencies for various equations of state 
of compact stars. These dashed lines represent the r-mode in- 
stability limits, and dotted lines the Keplerian limits. The cor- 
responding thick solid lines are the genuine upper frequencies. 
The light yellow district refers to the mass prediction of the X- 
ray pulsar in SAX J1808.4-3658 at 99.7% confidence lever ana- 
lyzed by Leahy et al. (2007), and the green dash-dot-dot horizon- 
tal line the corresponding 2.49 ms rotation frequency discovered 
by Wijnands and van der Klis (1998). 
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